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Abstract 

We analyse scalar flavour-changing neutral currents of down-type quarks in models with 
two Higgs doublets, coupled separately to up- and down-type quarks, in the limit where 
the ratio of the two expectation values (tan j3 = Vu/vd) is large. We clarify the origin of this 
phenomenon, both in AF = 1 and AF = 2 processes, analysing differences and analogies 
between supersymmetric and non-supersymmetric models. We conflrm previous flndings 
of a sizeable enhancement at large tan f3 of speciflc AF = 1 and AF = 2 amplitudes in the 
MSSM and, in these cases, we discuss how large-tan /? corrections can be controlled beyond 
lowest order. Finally, we emphasize the unique role of the rare processes Bs^d t^t~ 
and Bs^d — ^ /i'*'//" in probing this scenario. 



1 Introduction 



Processes mediated by flavour- changing neutral-current (FCNC) amplitudes are extremely 
useful to deeply investigate the dynamics of quark-flavour mixing: the strong suppression 
of these transitions occurring within the standard model (SM), due to the absence of 
tree-level contributions and the hierarchy of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, naturally enhance their sensitivity to possible non-standard phenomena. This 
is even more true in the case of down-type FCNC amplitudes mediated by scalar (and 
pseudoscalar) currents that, within the SM, are additionally suppressed by the smallness 
of down-type Yukawa couplings. 

The suppression of down-type FCNC scalar operators becomes less effective in models 
with an extended Higgs sector and, particularly, in the popular two Higgs doublets model 
(2HDM) of type II, where two SU{2)l scalar doublets are coupled separately to up- and 
down-type quarks (see e.g. Ref. In this case FCNC amplitudes are still absent at the 
tree level. However, contrary to the SM, it is possible to accommodate large down-type 
Yukawa couplings, provided the ratio Vu/vd = tan/?, where w„ {vd) denotes the vacuum 
expectation value of the doublet coupled to up (down)-type quarks, is large. Then, for 
instance, AF = 1 operators like bRSiJiR^L, suppressed by two small Yukawa couplings, 
could receive a tan^ /3 enhancement with respect to the SM case. 

Even though the tan^ f3 enhancement could be as large as 10^ in the simple 2IIDM, 
scalar FCNC amplitudes barely reach the level of vector-type SM amplitudes, and this 
happens only in very few cases. On the other hand, the 2HDM of type II is particularly 
interesting being the Higgs sector of the Minimal Supersymmetric SM (MSSM) (see e.g. 
Refs. 0). As first noted by Babu and Kolda 0, and later confirmed in |^, the tan/3 
enhancement of scalar FCNC amplitudes can be much more effective in the MSSM than 
in the non-supersymmetric case. In particular, the one-loop coefficient of the operator 
bRSLpR/iL has been found to scale like tan^/?, leading to possible large non-standard 
effects in Bg More recently, it has also been shown that the coefficient of the 

AF = 2 operator hRSLbiSR receives, at the two-loop level, a contribution scaling as tan^/3 
that could have a relevant impact on Bg-Bg mixing 0. 

The purpose of this paper is to clarify the nature of this tan/? enhancements. We 
will analyse the large-tan/? behaviour of all relevant down-type AF = 1 and AF = 2 
amplitudes, both in supersymmetric and non-supersymmetric cases. As already pointed 
out in 1^, the large-tan f3 enhancement of scalar FCNC amplitudes is intimately related to 
the appearance, at the one-loop level, of an effective coupling between and down-type 
quarks 0. As we shall show, this is a necessary consequence of any 2HDM and indeed 
it is realized in a very similar way in the supersymmetric and in the non-supersymmetric 
case. The difference between the two scenarios arises only under a specific conspiracy of 
the soft-breaking terms, which could decouple the supersymmetric Hu-Hd mixing from 
the ordinary H^-Hd coupling, necessarily suppressed in the large-tan/? limit. Analysing 
both scenarios, we shall clarify the origin of the different tan/3 factors, showing how to 
control the large-tan f3 terms beyond lowest order. 

The paper is organized as follows. Section 2 is devoted to AF = 1 amplitudes; there 
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we shall first analyse the generic structure of the effective di — dj — H^ vertex, then we shall 
discuss the large-tan (5 behaviour of di dj£^i~ amplitudes. In section 3 we analyse -6^,^- 
Bg^d mixing, discussing the generic structure of both reducible and one-particle irreducible 
contributions and deriving phenomenological bounds from AM^ rf. Section 4 contains a 
phenomenological analysis of the rare decays Bg ^ — > r+r^ and Bs ^. fi~^fi~ . The results 
are summarized in the conclusions. 



2 AF = 1 scalar currents 

2.1 The effective down- type Yukawa interaction 

In a 2IIDM of type II (including the MSSM) the tree-level Yukawa interaction is defined 
as 

C'y = dnYdQLHd + URYuQ lHu + h.c, (1) 

where Y^^d are 3x3 matrices in flavour space and Hu d denote the two Higgs doublets. 
This Lagrangian is invariant under a global U{1) symmetry, which we shall call U{l)d and 
under which da and Hd have opposite charge and all the other fields are neutral. If this 
symmetry were exact, the coupling of Hu to down-type quarks would be forbidden also 
at the quantum level. This symmetry, however, is naturally broken by terms appearing 
in the Higgs potential and, if tan (3 is large, this has a substantial impact on the effective 
Yukawa interaction of down-type quarks. Under the assumption that Ly is the only 
source of flavour mixing and that the mo del- dependent Higgs self- couplings have the 
MSSM structure |^, the one-loop effective down-type Yukawa Lagrangian (in both the 
supersymmetric and non-supersymmetric cases) can be written as 10) S 0: 



C 



eS _-] 



dnYd 



Hd + (eo + eyY^Yu) H:] Ql + h.c. , (2) 

(9(l/167r^)] proportional to the 



where eo,y denote appropriate loop functions [eo,y ~ 
C/(l)rf-breaking terms. 

In order to diagonalize the mass terms generated by C^, it is convenient to rotate the 
quark fields in the basis where Yd is diagonal [{Yd)ij = yf^ij]. In this basis we can write 



C 



eff 



Vddj^yt (l + eotan/5)5,, + ey tan/3 l^,";(|/i;)VO d^ + h.c 



(3) 



where t;^(l + tan^ /?) = {2^/2GF)~^ ^ (174 GeV)^ and denotes the tree-level CKM 
matrix, i.e. the CKM matrix in the limit ey = 0. Neglecting the small terms due to y'^2J 
defining Xjj^ = V^^V^/. (for j ^ k) and yt = y^, we can rewrite Eq. (|^) as 



•^d— mass 



Vdd'jiyf 



eyy^ tan /5 



-A 



4 + h.c. 



(4) 



where 



l + tan/? eo + eyZ/'|^.3 



0|2 



1 + tan/3 (eo + eyyf6, 



i3 



(5) 
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Because of the hierarchy of the CKM matrix, A* • <^ 1 and we can diagonahze Eq. 



perturbatively in A 
given by 



The rotation that diagonahzes Eq. 



to the first order in A*,- is 



di 



Sjk + ^YVt tan/5-^^ 



yjVj+yiyt ,t 



d\2 



Sjk + eyyt tan (3 



yht + ytyf 



(6) 
(7) 



At this order the eigenvalues are not shifted, and thus the leading tan /3 corrections to the 
eigenvalues of the down- type Yukawa matrix are simply described by Eq. (^) . On the 
other hand, the rotation (j^) modifies the structure of the CKM matrix, but only the V^j 
(j 7^ 3) elements receive 0{1) corrections (i.e. corrections not suppressed by the CKM 
hierarchy) 

1 + eotan/? 



1 + tan/?(eo + ey?//' 



As first noted by Babu and Kolda 0, if ey 7^ the diagonalization of Eq. (H) nec- 
essarily induces a FCNC coupling between quarks and neutral Higgs bosons. This can 
easily be understood by looking at the neutral component of Eq. (H): if ey 7^ and 
is not aligned to Yd, we have two independent flavour structures (not simultaneously di- 
agonalizable) that are weighted differently for physical Higgses and mass terms. Indeed 



performing explicitly the rotations 
terms, we find 



and (0) and neglecting the subleading 0{yj /y^^j) 



''-FCNC {k=ij) 



eyyt tan P 
1 + eo tan P 



tan P 



TjO* TjO 



+ h.c. 



(9) 



where, expressing the off-diagonal coupling in terms of the physical CKM matrix. 



^jk 



kj 



VkWs, 



l+tan/3(£o+£y»/g) 



1+eo tan /3 



l+ep tan /3 



1+tan l3{eo+eYyj) 



(k ^ 3, 3) 
(k = 3, 3) 
(k ^ 3, J = 3) 



(10) 



Equation (^ generalizes to any d'' — > d^ transition the result obtained by Babu and Kolda 
for the h s case. 

One of the most interesting aspects of Eq. (^ is the rapid growth with tan (3 of the 
H^dfi^di coupling. Expressing in terms of down-type quark masses, the FCNC coupling 
grows almost quadratically in tan (3, contrary to the approximate linear growth of the 
diagonal H^d^d\ term. This is the origin of the approximate tan^ (3 and tan^ (3 behaviour 
of the AF = 1 and AF = 2 amplitudes discussed in 0] and 0, respectively. 

As expected, the combination of neutral Higgs fields that appears in Eq. (^ has a 
vanishing vacuum expectation value and does not include the Goldstone component; it 
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Figure 1: Penguin diagrams generating a non- vanishing ey: a) non-supersymmetric 
Hd mixing; b) chargino mixing. 



can therefore be expressed in terms of the three physical neutral Higgs states. Assuming 
the Higgs potential to be CP-invariant and employing the notation of H we find 



^YVt 



cos [1 + eo tan /5] 



X 



cos(a - l3)h" + sin(a - - iA" (f^d^^ + h.c 



(11) 



2.2 Explicit estimates of ey 

In the context of a non-supersymmetric 2HDM with Higgs self-couplings fixed as in the 
MSSM case, the only source of U{l)d breaking is the bilinear operator HuH^- If the Higgs 
potential is CP-invariant, the coupling of this operator, conventionally denoted by yuP, 
is related to the mass of the neutral pseudoscalar Higgs by the (tree-level) relation 



/iP = sin 2/5/2 



(12) 



Thus U{l)d breaking is unavoidable if we require a non- vanishing mass for v4°. Note, 
however, that in this case the P(l)d breaking is parametrically suppressed in the large- 
tan /3 limit. 

The non-supersymmetric P^-P^ mixing induces a non- vanishing ey via the mechanism 
shown in Fig. |l]a. The computation of this effect at large tan (3 can be performed by setting 
(7 = 0, i.e. switching off gauge interactions: keeping g ^ Q would have complicated the 
calculation, owing to the gauge- dependent mixing of charged Higgs fields and W bosons, 
without affecting the final result. In this limit we can identify with the massless 
Goldstone boson and Pj" with the massive charged Higgs, and the only diagram to be 
computed is the one in Fig. |l|a, which leads to 



^2HDM 
-Y 



167r2 M^, 



log {m'jMl+)+0{m'jM^ 



H+ 



(13) 
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In the supersymmetric case an additional source of U{l)d breaking is provided by the 
fiHuHd term. The latter contribute to ey via the diagram of Fig. |I]b, leading to 0|: 

tL 

where, as usual, A denotes the coupling of the soft-breaking trilinear term (we assume 
both /i and A to be real), x^l = /^^/^f^ ^RL = ^jj^'f^- '^^^ expression of 
f{x,y) can be found in the appendix and the normalization is such that /(1,1) = 1/2. In 
the supersymmetric case one finds also a non- vanishing eo, dominated by the contribution 
of gluino penguins 0. 

As can be noted, the explicit expressions of efP^^ and ey^^^ are rather similar: ey ~ 
l/(167r^) times an adimensional coupling parametrizing the U{l)d breaking. In the non- 
supersjTumetric case the U{l)d breaking is strongly constrained by Eq. (jl^), which forces 
ey to be suppressed as 1/ tan/5 in the large-tan /3 limit. Thus the potential tan^/5 growth 
of the H^pfd^j^ coupling cannot be realized in the simple 2HDM model. On the contrary, 
in the supersymmetric case we are allowed to consider a scenario where nA/M~ = 0{1) 
also at large tan f3. Note, however, that this scenario implies a sizeable hierarchy among 
soft-breaking terms: if Mj^ > Ma we need A/B > tan/?, while if A/B < 1 we need 
M£ /Ml < tan/3. 

2.3 di ^ dj£~^£~ transitions 

We are now ready to discuss the effect of scalar-current amplitudes in di — > dji~^i~ 
transitions. The effective Hamiltonian describing these processes, including scalar-current 
operators, can be written as 

nfp^, = nil + CsOs + CpOp + c'sO's + c'pO'p + h.c, (is) 

where T^g^ denotes the SM basis of AF = 1 operators (see e.g. 0) and 

Os = IpdiU , Op = IpdiHi , 

O's = d[d^U , O'p = d[dih5^ . (16) 

One-loop contributions to b s,d transitions in the non-supersymmetric 2HDM 
have been discussed in detail in Refs. [^, |, |Tl|. This analysis requires, apparently, the 



evaluation of a large number of box and penguin diagrams. However, the calculation can 
be strongly simplified by working in the limit g = 0, which has proved to be the most 
convenient framework to discuss leading tan f3 contributions to scalar FCNC amplitudes. 
In this limit it is easy to realize that box diagrams can be neglected, as long as we ignore 
terms suppressed by two lepton Yukawa couplings. 

Penguin (and self-energy) diagrams can be evaluated by means of the effective Yukawa 
interaction discussed before, as shown in Fig. The limitation of the latter is that only 



5 



the leading dimension-four operator H*dfiQL has been considered. Higher- dimensional 
operators with more powers of Hu, which contribute to the mass diagonalization when 
Hu acquires a v.e.v., have been neglected. The coefficients of these higher-dimensional 
operators are necessarily suppressed by the inverse power of the heavy scale of the the- 
ory, namely the charged Higgs mass in the non-supersymmetric 2HDM. Therefore, using 
'^FCNC' control only the leading term in an expansion in powers of Vu/Mfj+. This 



leads to the following initial conditions for the Wilson coefficients in (|T5|): 



Cs — Cp 



— —Op 



327r2 

yfyevt^i 

327r2 



Ml. 



Ml. 



log 



log 



ml 

' Ml, - 



+ 



+ 



Ml., 



mt 



Ml., 



(17) 



where subleading 0{1/ tan/?) terms have been neglected. These results are compatible 
with those in Refs. ^ and, even though 0{m'f / Ml+) corrections are missing, 
Eqs. (0) have the advantage of incorporating higher-order tan /3 terms, hidden in yf and 
Xlj [see Eqs. (^ and (p!0|)]. As long as we neglect terms suppressed by additional down- 
type Yukawa couplings and/or additional off-diagonal CKM factors, these are the only 
source of tan /3-enhanced corrections. 

It is worthwhile to stress that the subleading 0{m^/Ml+) terms in Eqs. ([1^) can also 
be computed in a rather simple way. To this purpose we only need to perform a complete 
diagonalization of the off-diagonal mass term generated at one loop, which is equivalent 
to computing the diagram in Fig. Qb (in the gaugeless hmit). As a result, we find that 
the term within square brackets in Eqs. (|1^) should be replaced by 



Ml. 



MI+ - mf 




(18) 



in agreement with the results of Refs. P, P|, P|. 

Full one-loop calculations in the supersymmetric case have been performed in Ref. 
1^, complementing the result of Ref. obtained by means of the effective Yukawa cou- 
plings. As expected, the latter provides an excellent approximation to the full result in 
the limit of heavy squark masses. The supersymmetric contribution to Cs = Cp obtained 
by means of C^^q can be written as 



Cc 



yfyeytKi 



327r2 



5M| [1 + eo tan /?] 



md^mimf Xj, tan^ /3 



nA tan (3 



47r2 1 -I- tan/3(eo + eyyt^: 



i3) 



Mi Ml[l + eotan/?] 



(19) 

f{x^,L,XRL) ■ (20) 



The comparison between ([T9|) 
Some comments are in order: 



and (pOD illustrates the origin of the various tan (3 factors. 
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Figure 2: Feynman diagrams relevant to the c?^ (I-'l^l^r amplitude: a) leading contri- 
bution computed by means of the effective Yukawa interaction; b) one-loop contribution 
to the mass mixing in the 2HDM. 



1. The third tan (3 factor in the numerator is not directly related to the down-type 
Yukawa couplings. As anticipated, with a different choice of parameters it can be 
eliminated in favour of the ratio A/B. 

2. Higher-order corrections enhanced by tan/?, shown explicitly in Eq. (^), cannot 
be neglected (similarly to the 6 — ^ 57 case For tan (3 > 30 these higher-order 
terms are numerically more important than the remaining one-loop corrections not 
described by C'^^q. By means of Eq. ( pOD we take into account all terms of the type 
{Gp^hf^d X (a/7r)"'+-'^(tan/3)"'"'"^, where n > and a denotes either as or y^/^Air). 

3. The tan /3-enhanced corrections to the tRhLH^ vertex, which play an important 
role in 6 — > S7 although formally sub leading [0, are not relevant here owing to 
the strong overall suppression of the charged-Higgs contribution with respect to the 
chargino one. 

We conclude this section with a generic comparison between scalar FCNC amplitudes 
at large tan/5, both in the 2HDM and in the MSSM, and ordinary SM vector-type am- 
phtudes. On general grounds, the former are suppressed with respect to the latter by a 
factor 

where n ~ 2 in the non-supersymmetric case and n < 3 in the MSSM. An exception to 
this rule is provided by P — > decays, discussed in detail in section 4; here the helicity 
suppression of the SM matrix element leads to replacing the factor m^^m^ in (^) by m^,. 
Taking into account the approximate numerical relations 

and !I^,!^^of^') (22) 
mt rrit \tan^ p I rrit rrit \tanp J 
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that hold for 30 < tan/5 ^ 50, the general picture can be summarized as follows. 

b {s,d)T~^T~: non-standard scalar contributions could have a sizeable impact already 
at the level of inclusive transitions. The magnitude of the scalar amplitude could be 
substantially larger than the SM one in the MSSM, whereas it is typically smaller 
than the SM one in the non-supersymmetric case. 

b {s, d)fi^fi^: the impact of scalar currents is almost negligible in inclusive transitions 
(as exphcitly shown in Ref. 0), but large effects are still possible in the exclusive 
dilepton decays, especially within the MSSM. Indeed the most stringent constraint 
on AF = 1 scalar currents -under the hypothesis of minimal flavour violation- is 
at present derived from the experimental bound on B{Bs /^"""/i") (see section 4). 

s dfi~^ij,~: scalar current contributions could reach at most a few percent of the SM 



short-distance amplitude in Kl fJ^'^fJ. ■ Given the theoretical uncertainties affect- 



ing long-distance contributions to this mode [0, these non-standard effects are not 
detectable. We stress, however, that this conclusion holds only under the hypothesis 
of minimal flavour violation. If the CKM matrix is not the only source of flavour 
mixing, the eyfl parameters are not flavour- diagonal and the H^'c^j^d^ coupling is 
not necessarily proportional to Xjj. Within this more general framework, it is then 
possible to overcome the strong suppression of s — >• dfi'^fi~ amplitudes and to 
generate C(l) effects also in K decays. A detailed analysis of this scenario is beyond 
the purpose of this paper, where we shall restrict our attention to the hypothesis of 
minimal flavour violation and, henceforth, to B physics. 



3 Bd,s-Bd,s mixing 
3.1 Generalities 

Within new-physics scenarios where flavour mixing is governed only by the CKM matrix, 
the effective Hamiltonian relevant to mixing can be written as 



ttAB=2 '^F f./r2 



VtlVtidJ J2c^i^^)Q^ + h.c., (23) 



where the full basis of dimension-six operators can be found, for example in Ref. [|I3|. In 
addition to the SM operator 

Q"'"''' = h^f^qLh^qL , (24) 
those relevant at large tan /? are given by 

Q^^^ = InqLlnqL , 
Q^^"^ = bnqLhqR, 

qVRR = Inn^qRhRl^qR, (25) 
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Figure 3: Double-penguin diagram contributing to C^'"^; the black boxes denote the 
effective FCNC Yukawa interaction of Eq. (p!l|). 

where q = d,s. 

All the Wilson coefficients of the four operators in Eqs. (|2^ ) and (^) receives one- 
particle irreducible one-loop contributions from box diagrams. However, as pointed out 
in Ref . ^ , in the case of Q^^^ a sizeable contribution at large tan (3 is also provided by 
the reducible two-loop diagram in Fig. |^. 

The Wilson coefficients of the three operators in Eqs. (^) are all suppressed by small 
down-type masses, corresponding to the right-handed fields, but are possibly enhanced by 
appropriate tan (3 factors, as we shall discuss in the following. The smallness of implies 
that only Q^^^ plays a significant role in AM^. On the other hand, both Q^^^ and Q^^^ 
are potentially relevant to AMg. Despite a tan^ f3 enhancement, the contribution of Q^^'^ 
is practically negligible for both AM^ and AMg. 

Following the notation of Ref. the Bd^s^Bd^s mass difference can be written as 

AMd,s = |Vi(.,.)r , (26) 

where {Bb^^F"^ ) parametrize the hadronic matrix elements, to be determined through 
non-perturbative methods, and rjE is the QCD renormalization-group factor of the SM 
operator, given by 775 = 0.55 at next-to-leading order (NLO). The real functions -F/^'*, 
encoding short- distance contributions, are conveniently decomposed as 

Ftt' = S,{xtw)[l + hA , (27) 

where Xtw = '^V^w^ Sq^Xiw) accounts for the SM box (see appendix) and all non- 
standard effects are included in fd^s- Expressing fd,s in terms of the (non-standard) 
contributions to the Wilson coefficients of the four-fermion operators in Eqs. ( p^ and 
we can write 

C^^^(/ii) + C^^^(/it) + 4P^^^C^^^(/it) + 4P^^^C^^^(/it)l ' (28) 



45'o(xtvy) 
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where fit = 0{mt) and the exphcit expressions of the P\ taking into account renormahza- 
tion-group QCD corrections |13[ and matrix elements of the scalar operators -normalized 



to the SM one- can be found in Ref . ^ . Assuming 5-parameters of the scalar operators 
equal to one, setting mft(/ife) = 4.2 GeV, as{Mz) = 0.118 and Bb^s — 1-3 i 0.2 we find 



^pSLR ^ 3 5 ^ Q 5 ^ ^pSLL ^ i ^ Q 3 _ (-29) 

As can be noted, QCD effects enhance the contribution of scalar operators, especially 
qSlr^ thus even small new physics contributions to their Wilson coefficients may be 
relevant to the phenomenological analysis. 

3.2 = 2 box diagrams 

Computing explicitly non-supersymmetric box diagrams with the exchange of and if^ 
or two H"^ (Goldstone bosons included), in the large-tan/? limit, we obtain the following 
initial conditions for the Wilson coefficients 

VLL 2 

tan p 

= "^'Zhjf^ [LsixtH,XtH,l) -2L,{xtH,0,l) + Ls{0,0,l)] , 

^H^^ = JJ^ [LliXtW,XtwA) + Ll{XtH,XtHA) -'2Li{xtH,XtH,XwH)] , 

riSLR SmbTJid^stan^ (3 r . . \ , or f n ^ 

— 772 [—J^3[XtH,XtH,XwH)-r'^-l^3[XtH,^,XwH) 

—L^{0,0, Xwh) + XHwLi{XtH, XtH, xwh)] , (30) 

where, as usual, Xab = M^/M^ and the loop functions Li{x,y,z) are defined in the 
appendix. These results are in agreement with those recently reported in Ref. 0. 

The charged-Higgs contribution to C^^^ is directly suppressed by two inverse powers 
of tan/? and turns out to be completely negligible. The strong tan/3 enhancement of 
(jVRR jg jjiQj^g than compensated by the down-type quark-mass terms: employing the 
power counting in Eq. p^), Cjj^^ is effectively of C(l/ tan^/3) at most, and indeed it is 
numerically irrelevant. Similarly, also Cf/"^ is negligible being effectively of 0{1/ tan^/3). 
The potentially largest contribution -still suppressed with respect to the SM one- is 
provided by Cf/"^, which in the Bg-Bg case is effectively suppressed by only one inverse 
power of tan/5. For Mh ~ Cf/"^ can induce at most a 10% correction to F^^. 

^^2HDM-box ^ P^^^^Cl^^ ^ m,m tan^/3 r_^^^ _ ^^^^ _ ^ ^^^^^ _ 



2 

W 



^In principle the QCD correction factor of non-standard vector operators is not exactly 773; however, 
this difference can be ignored at the level of accuracy we are interested in. Similarly, we will neglect all 
finite QCD corrections to the initial conditions of the non-standard Wilson coefficients and corrections 
due to the running between /if and the new-physics scale. 
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Figure 4: Leading box-diagram contribution to C^^^ in the non-supersymmetric 2HDM. 



-0.10 - 0.07(xtH - 1) + 0{xtH - If 




(31) 



[mf,(/it) 3 GeV, ms(/ij) ^ 60 MeV]. We can therefore conclude that non-supersymmetric 
2HDM box contributions at large tan/5 do not induce appreciable effects in B^-B^ mixing. 

An important point to note is the following: whereas each down-type mass term is 
compensated by a corresponding tan (3 factor, as naively expected, in Cf/"^ and C^^^, 
this does not occur in Cjj^^ and Cf/"^. In the case of C^^^, the reason of the additional 
l/tan2/5 suppression is its vanishing in the absence of U{l)d breaking. Similarly to the 
case of Ai? = 1 amplitudes, this fact can easily be understood in the gaugeless limit. 
As shown in Fig. |: the two Hu-H^ mixing terms compensate the two tan (3 factors 
of the Yukawa couplings and, as a result, Cf/"^ has no explicit tan/3 dependence. As 
shown in section ^ U{l)d breaking is not necessarily suppressed at large tan/3 in the 
supersymmetric case, when the Hu-H^ mixing is replaced by the H^-Hd one. We can 
therefore expect a tan^ /^-enhanced chargino-squark contribution to C^^^. 

The full one-loop contributions to C^^^ and C^^^ generated by chargino-squark dia- 
grams in the MSSM can be written as 



(jVLL 

X 



C. 



SLL 



A[fiu,u)-2fiu,i) + fii,i)] 
4[g{i,i)-2g{u,i)+g{i,i)] , 



(32) 
(33) 



where 



fiu,i) 
g{u,i) 



i,j,h,k=l,2 



4 ^wxj irr- ^iuh Z^tk ^4 i^h 

i,j,h,k=l,2 ^"X] 



Xj' "^"hXi' -^XiXj/ ' 



(34) 
(35) 



and Yju^ and -Zj^^ are defined as in Ref. P^.p| 

The chargino-squark contribution to the Wilson coefficient of the standard vector 
operator (C^^^), which includes the supersymmetrization of the ordinary W-hox diagram, 

^ The superscript in Zj'-^ refers to the down-type Yukawa coupling involved. 
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has no special features at large tan/?: it is not suppressed or enhanced. For this reason, 
and also because C^^^ has been extensively discussed in the literature (see, in particular, 
Refs. [|l4l |15|1), we will not discuss it further. 



Contrary to C^^^, the chargino-squark contribution to the scalar operator Q^^^ has 
not been discussed before and, as anticipated, it has an interesting behaviour at large 
tan/3. In order to provide a compact analytical expression and a numerical evaluation of 
C^^^, we employ the following simplifying assumptions: 

• we neglect up and charm Yukawa couplings and any non- degeneracy among the first 
two generations of squarks: Mg. = M for q = u,c and i = 1,2 [as already implicitly 
assumed in Eqs. (p^)-(|35|)]; 



• we neglect off-diagonal terms in the chargino mass matrix (or we assume Mw <C 

max{M2,/i}), so that M^^ ^ M2 and M^^ ~ A*! 

• we neglect the left-right mixing of the squarks, but in the stop sector; there we 
introduce a mixing angle {—7t/2 < 9i < -k /2) satisfying the relation 

2mt{A — iicotB) 2mtA 
sin = T79 — ..9 ~ „ \ , (36) 

where the ii mass eigenstate (M^^ < M^J can be identified with the right-handed 
stop for small 9^ . 

In the limit where M^^ is degenerate with the first two generations of squarks, while Mj^ 
is kept lighter, we obtain the following rather simple expression: 

" 32G^^M^^/i^ <^(3^fi/.>3;f-/.) [giuj) = g{u,u) =0] , (37) 



where 



y2 



[x-yY + 0{x-yY , (38) 



and G2{y) can be found in the appendix. The expansion of G{x,y) around x ^ y lets us 
further simplify Eq. (^Tf), eliminating sin2^j in favour of A. We then obtain 

M^/i4 [1 + tan /3(eo + ey?//)] \MjJ 

where the Yukawa couplings have been expressed in terms of quark masses. Note the 
analogy of Eq. (ROT) to Eq. (p!^): in both cases the overall effect is ruled by the ratio 



fiA/M~ , which appears squared in the AB = 2 amplitude. However, since the effective 
FCNC interaction of Eq. (|^) does not enter the box amplitude, the tan j3 enhancement of 
the latter simply scale with the power of down-type Yukawa couplings. As we shall discuss 
later on, for large tan/?, small (/^ < 200 GeV) and fiA/M~ = 0{1), this supersymmetric 
scalar contribution has a small but non-negligible impact on AMj_ and AMg. 
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3.3 Double-penguin diagrams 



Double-penguin diagrams are formally a higher-order (two- loop) effect; however, the large- 
tan /3 growth of the effective Yukawa coupling make them numerically competing with 
one-loop contributions in the Bg-Bg case [^. 

The structure of Eq. ([11]) leads, in general, to a vanishing result when the Higgs 
fields are contracted to generate effective operators of the type bRqhbRqL or bL<lRbLqR 
{q = s,d). On the other hand, a non- vanishing result is obtained for effective operators of 
the type bnqLbiqR, or when the combination of Higgs fields explicitly shown in Eq. ([Til) 
is contracted with its Hermitian conjugate. As a result, only C^^^ is affected by double- 
penguin contributions. 

In general, the double-penguin contribution to C^^^ at large tan P can be written as 



C 

where ||6|1 



SLR _ ^T^^ybVs, 



1 + Co tan P 



1 + tan/?(eo + eyyt) 



cos^(q; — /3) sin^(a — /5) 



Mlo ^ M|o ^Mlo ^"mI" ^^^^ 

Within the non-supersymmetric case, the 1 / tan (3 suppression of e^H^^ leads to negligible 
effects for both AM a and AMg. On the contrary, the supersymmetric contribution to 
AMg can compete with the SM one. Expressing Yukawa couplings in terms of quark 
masses, we find 

fjSLR ^ GFrnbUismf tan^ pj^+p{x^L, xrl) ( jiA\ 

^ v^ttW^, [l + eotan/?]2[l + tan/5(eo + eyi/|)]2 l^Mi j ' ^ ' 

The tan^ j3 enhancement of this coefficient is partially compensated by the suppression 
and the additional l/lGvr^ (due to the two-loop order). As a result, it is of the same order 
as the chargino-stop box contribution to C^^^. Neglecting higher-order tan/3 terms one 
finds 

_ 47r^ mkMlh{x^L,xiiL) 

2 ^ 2 / s 2 

where the function h{x^L, xrl), normalized to h{l, 1) = 1, is plotted in Fig. ^. 





3.4 Numerical bounds from AM^^ in the MSSM 

In the case of AMd, the only relevant AB = 2 scalar operator is Q^^^. Neglecting possible 
new-physics effects in the standard vector-current operator, according to Eqs. (|37D and 
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Figure 5: h{Xfj,L, xrl) as a function of (= l/x^^) for different values of xrl. 
8D the large-tan P contribution to fd is given by 



fd 



Soixtw) Ml, 



M? 

tL 



Imposing that fd varies within the range allowed by the measurement of AMd, we can 
derive a set of bounds on the quantity fiA/M^^ for different values of stop masses (M^^ 
and Mf:^) and /i, as shown in Fig. ^ Although the experimental determination of AMd is 
rather precise, once the uncertainties on Bb^F^ and \Vtd\ are taken into account, fd can 
reach values close to 1, which we take as a reference figure. As can be noted from Fig. ^ 
the bounds on fiA/mf^ thus obtained are not severe. Even for small {fi ^ 100 GeV), 
light (m^^ ^ 200 GeV) and a large stop splitting, the upper bound on fiA/m^^ is 
above unity, condition that is naturally satisfied within our scenario and that still allows 
large effects in AB = 1 processes. 

Constraints even less strict on nA/M^^ can in general be set by using the present 
experimental lower bound for the Bg-Bg mixing. In this case also the double-penguin 
contribution to Q^^^ has to be taken into account. As shown in Eq. (^), for light fi and 
Ma ~ /i, the large tan/? contributions to C^^^ and C^^^ are comparable and tend to 
cancel out in fs, owing to the opposite sign between P^^^ and P^^^. Only if Ma is light 
and yU ^ Ma, or when the double-penguin contribution completely dominates, does one 
obtain a non-trivial bound. 
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Figure 6: as a function of fiA/Mf^ for /i = 100 GeV and different stop masses. The 
three sohd (dashed) curves from the left are obtained for {M£^(GeV), Mf^(GeV)} = 
{200, 500 (300)}, {300, 1000 (500)}, {500, 1000 (800)}. 



In the hmit where the double-penguin contribution dominates, higher-order tan/3- 
enhanced corrections are numerically relevant and cannot be ignored. In particular, terms 
proportional to |^ 

2as /iMg 



Co 



Ml 

dL 



(44) 



turn out to be very large. In order to obtain a simple expression, we can neglect the 
contribution of eyt/^ with respect to eo in the denominator of Eq. (|5), and set eo ~ 1/100, 
as obtained in the case of degenerate sparticles. We also set eo > 0, following the indication 
fi > from b ^ sj and {g — 2)^ |l^. We then obtain 



-0.4 



3(tan/?/50) 
1 + (tan/5/50) 



200 GeVy 
Ma ) 



M2 

tL 



(45) 



Since the present lower bound on AMg implies a lower bound on fs of about —0.5 §], 
we conclude that also in this scenario values of fiA/Mf^ < 1 are still allowed. Note that 
a much stronger bound would have been obtained without the inclusion of higher-order 
terms: for tan/5 ~ 50 these reduce the lowest-order result almost by a factor 5. 
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4 Phenomenology of Bs^d ~^ decays 
4.1 Effective Hamiltonian 

Within the SM there is only one dimension- six operator generating a non-vanishing con- 
tribution to Bq decays, namely 

Qio = hYqJl,l5i ■ (46) 
Being scale-invariant and completely dominated by short-distance dynamics, its Wilson 



coefficient is well approximated by the lowest-order result 



Lo - vr^ 1,T66G;V ) ■ 

The full NLO expression of Cf^ can be found in |p!8| . 

The additional scalar interactions that appear in the SM extensions under investigation 
are fully described by the operators in Eq. (|TBp. If the lepton Yukawa coupling is not 
negligible, as in the case of a t~^t~ final state, in principle we should take into account 
also the right-handed vector-current operator Q[q, obtained by Qio under the exchange 
{Qyb)L {q,b)ji. However, as can be expected from a naive extrapolation of the results 
in the previous section, the numerical impact of this operator is always negligible. Within 
the 2HDM one finds 

C,o - ^2 4M2,M^ ' ^^^^ 

where Bq{x) = [L3(x,0,l) — L3(0, 0, l)]/4. In the most favourable case, namely the 
Bg —>■ T^T" decay, leads at most to (9(1%) corrections to the SM result. A similar 
suppression is found also for the MSSM contribution to C[q. 

For completeness, we report here also the MSSM chargino-box contribution to Cs,p, 
which was not explicitly shown in section 0. This is given by 

Gs,P = 27r2 ^3i^3j X! ^l^Xn[^mth^m%^n!>fc^n4-^3(a;t;feXn5 ^^t^Xn' ^XmXn) 

h,k,m,n 

Tfl £ b ~ 

^ 7I~^^m4^mi>fc^'^'^fc^n4-^4(a;i>;,Xn5 3;t^Xn'^XmXn)] - ^ , (49) 

where we have assumed diagonal chargino-lepton-slepton couplings, and terms suppressed 
by rris^a/mh or by the mass difference between the ffist two generations of squarks have 
been neglected (under these hypotheses Eq. (^) is agreement with the results of Ref. p). 
Following the simplifying assumptions made for the chargino-squark contributions to -B"- 
B^ mixing, the above result becomes 



X _ ^x _ 9zMkv*V '^bmetan'' (3 r ^ 

- Gp - . 2 ^3i^3j 5 COS titL3[Xi^ 
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From this expression it is straightforward to check that, for large tan/5, box contributions 
cannot compete with the neutral Higgs-penguin coefficients in Eq. (pOD, as already antici- 
pated in section ^. We shall therefore neglect such terms in the following phenomenological 
discussion. 



4.2 Branching ratios 

Writing hadronic matrix elements of axial and pseudoscalar currents as 



q 1 



{Q\q-f,h\B,{p)) = -ifs. 



Ml 



(51) 



the most general expression of Bq — > 
malization of the effective operators: 



{rrib + mg) 

branching ratios reads, according to our nor 



fLMBjB.mj 



1 - 



67T 



4mj 



1/2 



1 - 



4m| 
Ml 



MUCs-C's) 



2mi{mb + ruq) 



+ 



MUCp-C'j 



2mi(mb + rrin 



10 



(52) 



where light-quark masses and Wilson coefficients are understood to be evaluated at a scale 
of C(mfe). Since the ratios {Cs,p — Csp)/me are independent of mi, to a first approxi- 
mation the relative weight of the various contributions is independent of m^, whereas the 
overall branching ratio scales like mj. This leads to a strong enhancement factor of the 
r+r^ modes that partially compensate their difficult detection. 

Within the SM, employing the full NLO expression of Cf^ 
of the two Bg modes are given by 



Tsj , the branching ratios 



B{B, 
B{B, 



T T 



SM 



SM 



3.1 X 10"^ j 
215 X B{Bs 



Bs 



ts\ 



mt[mt) 



3.12 



0.21 GeV 



0.04, 



1.6 ps/ \166 GeV^ 



SM 



(53) 
(54) 



whereas the corresponding B^ modes are both suppressed by an additional factor \ Vtd/Vt 



ts\ 



= (4.0 ± 0.8) X 10-2. Note that the B {B, t+t-)/B{Bs yU+Zi") ratio in Eq. (||) 
can be considered as a model-independent upper bound: if Cs,p oc this ratio can 
only be smaller than in the SM case. At present the experimental bound closest to SM 
expectations is given by 



B{B, 



< 2.6 X 10-^ (95% CL) [19] 



(55) 



Working in the limit where the supersymmetric contribution to p in Eq. ( PP] ) domi- 
nates over all other contributions, and neglecting terms suppressed by m^/mfc and m^/m^. 
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we find 



MSSM 2567r5 V M 




lv(l + eotan/3)[l + tan/3(eo + eyy2)]j • ^ > 



Setting eo ~ 1/100 and neglecting the ey?/^ term in the denominator [as done in Eq. (^S])], 
and employing the reference values of /b^, tb^, rrit and Vts as in Eq. (^31), we obtain the 
following compact phenomenological expression 

BiB ^ u^u~) ^ 3 X 10-^ /200 GeVy / /i^W^^ ' 

/i j^^g^^dxlO Q^l^^J4(, Ma J \ Ml J ' 

where r = tan/?/50. Given the weak constraints on fxA/M- discussed in section |], we 

cannot exclude a scenario where B {Bs —>■ fJ'~^fJ'~) is just below its present experimental 
bound, as already pointed out in Refs. |P, H- However, we stress that this can occur 
only for rather extreme values of the parameter space, in particular for tan/5 ~ 50, 
when higher-order tan /3-enhanced corrections are very large. For tan/3 < 30 it is already 
unlikely to find B {Bs above the 10~^ level and for tan/5 < 10 supersymmetric 

contributions do not exceed the level of the SM ones. Concerning higher-order tan/3- 
enhanced terms, we recall that these are fully taken into account by means of Eq. ([56|), 
once the appropriate expressions of eo and ey are used, whereas the approximate result 
in Eq. (H) is valid only for M3 ~ M,- ~ /i > 0). 

In the scenario where the present discrepancy between experimental data and SM 
predictions on the anomalous magnetic moment of the muon |^0|, ^ is due to super- 
symmetric effects, the possibility of a sizeable enhancement of B {Bg —>■ fi'^fi~) certainly 
becomes more likely. Indeed in both cases supersymmetric contributions grow with tan (5. 
As recently pointed out in Ref. [Q, the correlation between these two observables be- 
comes very strong in the constrained minimal supergravity scenario. Unfortunately, the 
situation is not so clear in the general framework of the MSSM. A useful tool to illustrate 
basic features of the supersymmetric contribution to = {g — 2)^j/2 is the expression 



15 X 10-10 ■ V 50 ; V / 



which provides a good approximation to the full one-loop result [|16| in the limit of almost 
degenerate higgsinos and electroweak gauginos (Mi ~ M2 ~ ;U ^ Mw), and allowing a 
moderate splitting between slepton and chargino masses (the normalization of the l.h.s 
is the present experimental error l2^). Comparing Eqs. (|57D and (|58D helps us draw the 



following general conclusions, whose validity goes beyond the approximations made to 
derive the two equations: 
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The mild tan j3 dependence of aj^^^^ allows significant supersymmetric contributions 



to the latter, even for tan (3 < 10. In this case the effects on B {Bg jji^ jj, ) would 
be undetectable, at least before the LHC. 



B{Bs strongly depends on A and Ma, which play almost no role 

in aj^^^^. As a result, B{Bs /x"*"/!")^^^^ can be made arbitrarily small even if 
qMssm saturates the experimental result and tan /? ~ 50. 



• A detection of B{Bs 1^^ at the 10"^ level (or above) would provide a strong 
support in favour of the supersymmetric explanation of the {g — 2)^ puzzle. This 
result could be interpreted in the MSSM framework only for tan/? > 30 and light 
Ma, a scenario where it is difficult to keep aj^^^^ below 10^^ unless chargino masses 
are unnaturally high. 

Concerning Bs^ —* t^t^ modes, the following relations hold in the limit where scalar- 
current contributions are dominant: 



B{Bs T+T-) _ 

B{Bd^T+T-) 



(59) 



Bs 



^ Bii ^B^ 


Vtd 


^tb, ml 


PbMI 


Vts 


^tb, ml 




B^B^-.^^^-^ .2 ,.5 ^ ..2 ^^ ..2, 6.8 ± 1.4 . (60) 

By means of Eq. (^), these allow us to set the following indirects limits: 

B {Bs r+r") < 4.3 x 10"^ , B {Ba ^ r+r") < 2 x 10^^ . (61) 

Unfortunately the present experimental information on the r+r^ modes is very poor: we 
are not aware of any progress with respect to the indirect bounds, at the per cent level, 
discussed in Ref . . The signature of these modes is certainly very difficult at hadron 
colliders, but significant improvements could be expected from B factories. Reaching a 
10~^ sensitivity on the B^ — > r+r" mode, 5-factory experiments could compete with the 
present search for Bg l^^ at Fermilab. 



5 Conclusions 

In this paper we have presented a general analysis of scalar FCNC amplitudes of down- 
type quarks in two classes of models with two Higgs doublets: the MSSM and the non- 
supersymmetric 2HDM of type II. We have shown that the tan/3 growth of these am- 
plitudes is not simply determined by the number of Yukawa couplings appearing in the 
amplitudes. A crucial role is also played by the breaking of the U{l)d symmetry that 
forbids the tree-level coupling of Hu to down-type quarks. Within the MSSM, the U{l)d 
breaking induced by the Hu-Hd mixing is not necessarily suppressed in the large-tan/? 
limit, resulting in a potential tan^ (3 growth of the effective d'^^d^Hd vertex generated at one 
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loop. This anomalous tan f3 behaviour is characteristic of the effective FCNC couphng of 
down-type quarks to neutral Higgses and it does not appear in irreducible AF = 2 ampli- 
tudes at the one-loop level. On the other hand, because of the 1/ tan/3 suppression of the 
ordinary Hu-H^ mixing, this phenomenon does not arise at all in the non-supersymmetric 
model. As explicitly shown, in both cases leading tan (3 contributions to scalar amplitudes 
can efficiently be computed in the gaugeless limit of the models, or considering only the 
Higgs sector of the theory. 

Despite the potential large enhancements (if tan/3 > rrit/mh), non-standard scalar 
FCNC amplitudes are difficult to be identified experimentally. The rare dilepton decays 
Bs,d — > offer an almost unique opportunity in this respect. As pointed out first in 

Ref. 0, the tan^/3 growth of the b — >■ s/i"'"/i~ scalar amplitude, within the MSSM, could 
lead to order-of-magnitude enhancements of the rate compared to SM expec- 

tations. We have explicitly checked that this statement remains true even when taking 
into account the existing constraints on AB = 2 scalar operators. The latter includes 
two types of effects: i) the reducible two-loop contribution proportional to tan^ /? 
discussed in Ref. 0; ii) an irreducible one-loop contribution proportional to y^tan^/?. 
The two effects (comparable in size in the Bg-Bg case) are typically smaller than the SM 
amplitude and weakly constrained at present. 

The order of magnitude enhancements of Bg^d rates can occur only for tan (3 > 

30, when higher-order tan /3-enhanced terms cannot be neglected. Following the approach 



of Ref. |I0[, we have discussed how to control these terms beyond lowest order. Despite 
a drastic reduction of the one-loop result due to higher-order tan /3-enhanced terms (up 
to a factor 5 for tan/3 ~ 50), the possibility that B {Bg — > fi^fi^) is just below its present 
experimental bound is still open. An evidence of B{Bg f^^fJ'") at the 10~^ level (or 
above) would provide a strong support in favour of a supersymmetric explanation of the 
{g — 2)^ puzzle |2^; however, we have shown that the opposite is not true in the general 
MSSM framework: a sizeable supersymmetric contribution to {g — 2)n does not necessarily 
imply observable consequences in B{Bg fi~^fi~). Finally, we have emphasized the 
interest of Bg^d t^t^ modes in probing non-standard scalar FCNC amplitudes and, 
particularly, the importance to search for B^ t^t^ at B factories. 
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A Loop functions 



Soix) 

F{x,y,z) 

Li{x,y,z) 
L2{x,y,z) 
L3{x,y,z) 
L4{x,y,z) 

G2{y) 



Ax — llx"^ + x^ 3 lo£ 



X 



4(1 -x) 



2(l-xY 



So 



mi 
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X \ogx 
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